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Dynamical phase transitions in long-range Hamiltonian systems and Tsallis 
distributions with a time-dependent index 
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We study dynamical phase transitions in systems with long-range interactions, using the Hamil- 
tonian Mean Field (HMF) model as a simple example. These systems generically undergo a violent 
relaxation to a quasi-stationary state (QSS) before relaxing towards Boltzmann equilibrium. In 
the coUisional regime, the out-of-equilibrium one-particle distribution function (DF) is a quasi- 
stationary solution of the Vlasov equation, slowly evolving in time due to finite A*' effects. For 
subcritical energy densities, we exhibit cases where the DF is well-fitted by a Tsallis g-distribution 
with an index q{t) slowly decreasing in time from q ~ 3 (semi-ellipse) to q = 1 (Boltzmann). When 
the index q{t) reaches an energy dependent critical value qcrit, the non- magnetized (homogeneous) 
phase becomes Vlasov unstable and a dynamical phase transition is triggered, leading to a mag- 
netized (inhomogeneous) state. While Tsallis distributions play an important role in our study, 
we explain this dynamical phase transition by using only conventional statistical mechanics. For 
supercritical energy densities, we report for the first time the existence of a magnetized QSS with a 
very long lifetime. 
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Systems with long-range interactions have a pecuUar 
dynamics and thermodynamics [l|. In these systems, the 
potential of interaction decreases at large distances like 
r~", with an exponent a smaller than the dimension of 
space. As a result, the energy is non-additive and this 
can lead to negative specific heats and inequivalence of 
statistical ensembles. The equilibrium properties of these 
systems are now well-understood and the main challenge 
concerns the description of their out-of-equilibrium prop- 
erties. The general methodology to study these systems 
can have impact in different domains of physics (see var- 
ious contributions in [l|). 

Long-range systems can spontaneously develop out- 
of-equilibrium coherent structures, called quasi station- 
ary states (QSS). The QSS are generically characterized 
by non-Boltzmannian distributions and their lifetime di- 
verges with the size N of the system. Galaxies in the 
universe and 2D vortices in geophysical and astrophysi- 
cal flows are examples of such QSS [j]. They arc also 
observed in non-neutral plasmas under a strong mag- 
netic field [31 , in the free electron laser [^ and in simple 
models like the Hamiltonian mean-field (HMF) model 
[a, la, 0, Hi B ll3l- They are therefore general features 
of systems with long-range interactions. The interest- 
ing issues about the nature of these QSSs concern: i) 
their rapid formation starting from a generic non equi- 
librium initial condition; ii) their structure and stability; 
iii) their slow evolution due to finite N effects; iv) their 
final transition to the Boltzmann equilibrium distribu- 
tion. In particular, point iv) requires to develop a theory 



of dynamical phase transition adapted to these systems. 
In this perspective, the HMF model [11| became a 
paradigmatic laboratory. It represents a system of N 
globally coupled rotators of unit mass described by the 
Hamiltonian: 



H= - 



1 ^ 
1=1 



1 ^ 



2N 



e,)], (1) 



i,J = l 



where Oi represents the orientation of the i-th rotator and 
Pi is the conjugated momentum. The evolution of the sys- 
tem may be followed through the magnetization, an or- 
der parameter defined as M = |M| = | ^ va.i\/N , where 
m.i = (cos6'i, sin^i). The HMF model has been exten- 
sively studied as a representative of a broad class of sys- 
tems with long-range interactions, like gravitational sys- 
tems (see references in [7| ) . The equilibrium solution 11 1 
reveals the existence of a second-order phase transition 
at the critical energy density Uc = 3/4: the Boltzmann 
equilibrium state is homogeneous (non-magnetized) for 
U > Uc and inhomogeneous (magnetized) for U < Uc- 

On the other hand, the interpretation of QSS is an is- 
sue of tense debate that has recently polarized researchers 
in two groups. The first group @ interprets the QSS as 
non-Boltzmannian equilibrium states explainable within 
Tsallis generalized thermodynamics [l2| when the ther- 
modynamic limit N — > +00 is taken before the infinite 
time limit t — > -|-oo. The justification advocated is that 
the evolution of the system is generally non-ergodic so 
that non-standard distributions can emerge. The other 
group p, [tI, [SJ, |9[ interprets the QSS in terms of the 



Lyndcn-Bell [13[ theory of violent relaxation that was ini- 
tially developed for coUisionless stellar systems and 2D 
vortices l^- For sufficiently short times, or for N —> +00, 
the evolution of the one-particle distribution function 
(DF) is governed by the Vlasov equation, a mean field 
equation that ignores correlations between particles [l4 1 . 
In this context, the QSS result from a process of phase 
mixing and violent relaxation on a coarse-grained scale. 
The statistical theory of Lynden-Bell [l3| predicts the 
most mixed state consistent with all the constraints of 
the coUisionless (Vlasov) dynamics. 

Numerical simulations have given contrasting results. 
For a water-bag initial condition (i.e., uniform distribu- 
tion of velocities and angles in an interval symmetric 
around 0) with energy density U ~ 0.69, it has been 
found that the Lynden-Bell theory predictions are correct 
[9| if the initial value of the magnetization Mq = M{t = 
0) is below a critical value (Mo)crit = 0.897 appearing 
in the theory JsL l9|, while this is no more true if Mq is 
close to 1 [5|, llO|. In particular, for Mq = 1, Campa et 
al. [lO| find that the QSS is well- fitted by a homogeneous 
{Mqss ~ 0) semi-elliptical distribution, which is a par- 
ticular Tsallis DF with g = 3 (as noted in [ij]), while the 
Lynden-Bell prediction is an inhomogeneous {Mqss ¥" 0) 
Boltzmann distribution corres pon ding to the non degen- 
erate limit of the theory [ij, ll5|. Stimulated by these 
results, one of us (8[ has proposed an approach in which 
the two points of view can be reconciled in the follow- 
ing manner: (i) if the system mixes well (ergodicity), 
the QSS is (close to) the "most mixed state" predicted 
by Lynden-Bell's theory; (ii) if the system does not mix 
well (non ergodicity), relaxation is incom,plete and non- 
standard DFs (including the important class of Tsallis 
distributions) can occur. 

The purpose of this Rapid Communication is to give 
numerical evidence that Tsallis distributions can also 
characterize the slow collisional evolution of the QSS (due 
to finite N effects) up to the relaxation to Boltzmann 
equilibrium. Our results support the proposal made in 
Ref. [14|, according to which the out-of-equilibrium DF, 
for a large class of initial conditions, could be fitted by a 
Tsallis g-distribution (see below) with an index q{t) de- 
pending on time. We will show that the dynamical phase 
transition from the non-magnetized QSS to the magne- 
tized Boltzmann distribution is triggered when the value 
q{t) is such that the DF becomes Vlasov unstable: there- 
fore, the final transition to the Boltzmann distribution 
is driven by a dynamical instability with respect to the 
Vlasov equation. 

In a non-magnetized QSS, the DF depends only on 
velocity and slowly evolves in time due to finite N effects. 
It turns out that the time-dependent DF is very well 
fitted by Tsallis distributions (polytropes): 

l/(9(t)-l) 



with a time dependent index q{t) > 1. The support of the 
function is the interval [—Vmaxit),Vmaxit)]] the normal- 
ization factor C(i) and Vmax(t) depend on time through 
q(t). The normalization condition /_ fdv = l/(27r) 
and the conservation of energy U — 1/2 = tt J_ fv^dv 

lead to the relations l/(27r) = Cv,naxV^^{^ + Vl? ^ 
l))/r(3/2+l/(g-l)) and C/-1/2 = 7rCi;3_(V^/2)r(l+ 
l/(q - l))/{(r(3/2 + l/(g - l))(3/2 + l/{q - 1))} from 
which we obtain the identity v1^^^{t) = 4(C/- l/2)(3/2-f 
\/{q{t) — 1)). On the other hand, using the stability 
criterion 1 + t^ J_^{f'{v)/v)dv > of [y, [l^ (see also 
[7|) we find that a spatially homogeneous g-distribution 
is Vlasov stable (dynamical stability) iff t/ > Ucrit{q) = 
l/2+(l/4)(q+l)/(3q-l) 0,0, [13. Forg = l (Boltz- 
mann) we recover t/crit(l) = Uc = 3/4, for q -^ +00 
(water-bag) we recover [/* = Ucrit{oo) = 7/12 and 
for g = 3 (semi-elliptical), we recover C/crit(3) = 5/8 
[lO|. From this criterion, DFs of the form ^ arc al- 
ways stable for U > Uc ^ 3/4 and always unstable for 
t/ < t/* = 7/12. For U^ < U < Uc, they are stable iff 
q > qcvu{U) with [H 



qcrit{U) 



4C/- 1 
12C/-7' 



(3) 



f{v,t)=C{t) 



1 



cW 



(2) 



Campa et al. [lO[ considered a water-bag initial condi- 
tion with Mq = at energy density U = 0.69. This initial 
DF is already a maximum Lynden-Bell entropy state so it 
does not experience phase mixing and violent relaxation 
14| . However, it slowly evolves due to finite N effects. 
In Ref. [10[ it was found that the system rapidly forms a 
velocity DF of semi-elliptical shape. As we have already 
indicated, this is a Tsallis DF of the form ^ with q = 'i. 
However, the slow evolution of the DF on longer times 
was not quantitatively addressed. Here, we have per- 
formed simulations at several energies, all of them start- 
ing from an unmagnetized state (uniform angle distribu- 
tion). For the subcritical energy densities U < Uc, the 
initial velocity distribution was a semi-elliptical, q ~ 3, 
function. The energy densities were all above the criti- 
cal value 5/8 corresponding to g = 3, and therefore the 
initial steady states were Vlasov stable. In these cases 
we found that during the QSS, while the magnetization 
remains zero (or, more precisely, oscillates about a value 
scaling to zero as \/Vn because of finite size effects), the 
velocity DF has a slow evolution. At a certain point the 
distribution is no more Vlasov stable; then the magneti- 
zation begins to increase relatively fast, and the systems 
heads towards Boltzmann equilibrium. During the slow 
process characterizing the QSS, we found that the veloc- 
ity DF remains always close to one of the q functions ([2]), 
so that it is possible to define a function of time q{t). In 
this way, the evolution of the velocity DF is described 
by the slow decay of a single quantity, the index q{t) of 
the polytrope, from the value g(0) = 3 to the energy 
dependent critical value qcrit{U). 



In Fig. 1, we show results corresponding to U — 0.69 
(similar results, not shown, were obtained for U = 0.65, 
U = 0.67, and U = 0.71). The left panel shows the evo- 
lution of the magnetization M(t), while the right panel 
shows the evolution of the function q{t). The runs corre- 
spond to A^ = 2^^. The critical value of q corresponding 
to U = 0.69 is qcr»t(0.69) = 1.375. This value is indi- 
cated, for reference, as a dashed line in the graph of q{t). 
At each time, the best value of q is determined by a least 
square procedure, minimizing with respect to q the sum 
of the squares of the differences between the numerical 
distribution and the g-distribution for the energy density 
U = 0.69. In this procedure, both the numerical and the 
theoretical distributions are given as histograms with a 
suitably chosen bin width; q is varied in steps of 0.001. 
We see that the destabilization of the DF, and the associ- 
ated process of magnetization of the system, is related to 
the approach of q{t) towards the critical value qcrit{U)- 
We found that this behavior is present for all N values 
that we have studied, although noise obviously increases 
if N is decreased. 
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FIG. 1: Time evolution of the magnetization M{t) (left panel) 
and of the index q{t) (right panel), for U — 0.69 and A'' = 2^"^. 
The dotted line in the right panel indicates the critical value 
gcrit(0.69) = 1.375. During the coUisional process, the DF is 
a quasi-stationary solution of the Vlasov equation of the form 
^ slowly evolving in time due to finite A'^ effects (see Figure 
2). When q{t) reaches qcrit, it becomes dynamically (Vlasov) 
unstable and the phase transition is triggered. 



In Fig. 2 we show snapshots of the DF during the 
"coUisional" relaxation for the same run as that shown 
in the previous figure. The full lines give the numerical 
distributions at a given time, and the dotted lines are the 
g-functions for the best fitted value of q{t). The fit by 
a Tsallis DF with a time dependent index is very good 
during all the out-of-equilibrium evolution. 

Finally, we consider the case of a supercritical energy 
density U > Uc = 3/4. For U = 0.8, Campa et al. 
[lOJ I started from a non-magnetized state with a semi- 
elliptical distribution {q = 3), which is Vlasov stable, 
and found that the system relaxes to the Boltzmann equi- 
librium distribution on a very slow timescale scaling like 
tcoiiiN) ~ e^. Since qcrit{U) < 1 for [/ > U^ while Tsah 
lis distributions with a compact support have q{t) > 1, 




FIG. 2: Velocity DF at four different times during the QSS 
of the run shown in Figure 1. The full lines are the numerical 
distributions, while the dotted lines are the g- functions for the 
best fitted value of q{t). a) t = 3 x lO'' and q{t) = 2.224; b) 
i = 6 X 10^ and q{t) = 1.987; c) t = 9 x 10^ and q{t) = 1.789; 
d) f = 1.2 x 10" and q{t) = 1.697. The absolute value of the 
DF depends on the chosen bin width. 



the spatially homogeneous phase remains always Vlasov 
stable during the coUisional evolution and no dynamical 
phase transition occurs. This explains 1J| the long ob- 
served relaxation time for U > Uc ss compared to the 
case U < Uc- In this case, the relaxation is due to a 
slow coUisional relaxation, not to a dynamical (Vlasov) 
instability. Considering a slightly different supercritical 
energy density, U = 0.85, for which qcrit = 0.75, the left 
panel of Fig. 3 shows the function q{t) that, starting 
from q{0) = 3, slowly approaches 1. During this relax- 
ation the DF is well approximated, as in the subcritical 
case shown in Fig. 2, by the g-function with q = q(t). 
We have also studied another type of initial condition, at 
the supercritical energy density U = 0.85, starting from 
a non-magnetized state with a Tsallis velocity DF cor- 
responding to g = 0.5, which is below the critical value 
0.75. This initial DF is Vlasov unstable, linear insta- 
bility followed by phase mixing and violent relaxation 
quickly drive the system towards a magnetized out-of- 
equilibrium QSS. Although the system should eventu- 
ally relax to a non-nragnetized Boltzmann distribution, 
the right panel of Fig. 3, corresponding to a run with 
N = 2^^, shows that the duration of this magnetized 
QSS is extraordinary long (probably exponential in TV, 
as estimated by runs at different N values) . For compar- 
ison, we note that for TV = 2" the QSS at U ^ 0.69 lasts 
for a time of about 5 x 10^ [l0| . 



In conclusion, we have shown that TsaUis q- 
distributions (polytropes) can be useful to describe cer- 
tain features of the complex HMF dynamics. First, they 
can provide good fits of the QSS in certain cases of in- 
complete relaxation 10|,llJ]. However, this is not general 



14j and other DF that are stable stationary solutions 
of the Vlasov equation can be observed as well (the sys- 
tem can also have a non-stationary oscillatory behaviour 
[1911). More interestingly, like in the gravitational simu- 
lations of Taruya and Sakagami [20| , we have found that 
the out-of-equilibrium evolution of the DF in the coUi- 
sional regime can be fitted by Tsallis distributions with 
a time dependent index q{t). This may be a general fea- 
ture of the dynamics of long-range systems. However, 
while Taruya and Sakagami [20| interpret the dynamical 
phase transition as a generalized thermodynamical in- 
stability (in Tsallis sense), we interpret it as a d yna mical 
instability with respect to the Vlasov equation [ij, |21[ . 




sional relaxation time tcoii {N) is larger than Nto (where 
t]j is a typical dynamical time of the system) , but its pre- 
cise scaling with N cannot be deduced from the present 
kinetic theory. Different scalings have been found nu- 
merically depending on the initial condition: Nto [5||, 
N^-Hd U El and e^to [Si- Due to the considerable 
timescale separation between the dynamical time (of or- 
der tjy) and the collision time (larger than Ntn), it can 
be argued that, during the slow "coUisional" relaxation, 
the time dependent DF remains close to a steady solu- 
tion of the Vlasov equation that slowly evolves under 
the effect of collisions (finite A^ effects). In this Rapid 
Communication we have shown that, at least for a class 
of initial conditions, these slowly evolving QSS are well- 
fitted by Tsallis (7(i)-functions. However, as discussed 
in Sec. 5 of [2l| this is not necessarily the mark of a 
"generalized thermodynamics" [I4I . Indeed, if using only 
classical methods of kinetic theory (consistent with the 
Boltzmann entropy), we were able to write down the ki- 
netic equation describing the coUisional relaxation of the 
spatially homogeneous HMF model (the analogue of the 
orbit-averaged- Fokkcr-Planck equation in astrophysics), 
the DF ([2|) could form an approximate time dependent 
solution of this classical kinetic equation (without rela- 
tion to generalized thermodynamics). Unfortunately, for 
the HMF model, this kinetic equation is not known. 
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FIG. 3: Time evolution of the index q{t) (left panel) and of the 
magnetization M{t) (right panel), for U = 0.85 and A'' = 2^^. 
The left panel refers to a run with q{0) = 3, in which the 
system remains always demagnetized {q{t) > Qcrit = 0.75) 
and slowly approaches Boltzmann equilibrium [q = 1); the 
right panel refers to a run with q{Q) = 0.5 < qcrtt = 0.75, in 
which the QSS is magnetized and lasts for a very long time. 



The evolution of the system from the QSS to the Boltz- 
mann distribution, due to the development of correlations 
("collisions") between particles, could be understood by 
developing a kinetic theory of the Hamiltonian A^-body 
problem. As a matter of fact, starting from the Liouvillc 
equation and using an expansion of the solutions of the 
BBGKY-like hierarchy in powers of 1/N in a proper ther- 
modynamic limit A^ -^ +00, a general kinetic equation 
can be obtained at the order 0{\/N) [li|; the BBGKY 
hierarchy is closed, neglecting three body and higher cor- 
relations of order 1/7V^ or smaller. However, for one 
dimensional systems like the HMF model, this kinetic 
equation does not tend to the Boltzmann distribution 
(for spatially homogeneous systems the collision term, 
which reduces to the Landau or Lenard-Balescu opera- 
tor [ij, I24I, cancels out identically in \D). This implies 
that higher order correlations are needed to understand 
the convergence of the system to the Boltzmann distri- 
bution. Therefore, we can only conclude that the colli- 
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